Astronomy & Astrophysics manuscript no. shellsPNG 


© ESQ 2008 


February 5, 2008 





HST/ACS observations of shell galaxies: inner shells, shell colours 

and dust. * 

G.Sikkema 1 , D.Carter 2 R.F. Peletier 1 , M. Balcells 3 , C. del Burgo 4 , and E.A. Valentijn 1 

1 Kapteyn Astronomical Institute, University of Groningen, PO Box 800, 9700 AV Groningen, The Netherlands 

2 Astrophysics Research Institute, Liverpool John Moores University, 12 Quays House, Egerton Wharf, Birkenhead, CH41 1LD, 
United Kingdom 

3 Instituto de Astrofisica de Canarias, Via Lactea s/n, La Laguna, 38200, Spain 

4 School of Cosmic Physics, Dublin Institute for Advanced Studies, Dublin 2, Ireland 



ABSTRACT 

Context. Shells in Elliptical Galaxies are faint, sharp-edged features, believed to provide evidence for a merger event. Accurate pho- 
tometry at high spatial resolution is needed to learn on presence of inner shells, population properties of shells, and dust in shell 
galaxies. 

Aims. Learn more about the origin of shells and dust in early type galaxies. 

Methods. V-I colours of shells and underlying galaxies are derived, using HST Advanced Camera for Surveys (ACS) data. A galaxy 
model is made locally in wedges and subtracted to determine shell profiles and colours. We applied Voronoi binning to our data to 
get smoothed colour maps of the galaxies. Comparison with N-body simulations from the literature gives more insight to the origin 
of the shell features. Shell positions and dust characteristics are inferred from model galaxy subtracted images. 
Results. The ACS images reveal shells well within the effective radius in some galaxies (at 0.24 r c = 1.7 kpc in the case of NGC 5982). 
In some cases, strong nuclear dust patches prevent detection of inner shells. Most shells have colours which are similar to the under- 
lying galaxy. Some inner shells are redder than the galaxy. All six shell galaxies show out of dynamical equilibrium dust features, like 
lanes or patches, in their central regions. Our detection rate for dust in the shell ellipticals is greater than that found from HST archive 
data for a sample of normal early-type galaxies, at the 95% confidence level. 

Conclusions. The merger model describes better the shell distributions and morphologies than the interaction model. Red shell colours 
are most likely due to the presence of dust and/or older stellar populations. The high prevalence and out of dynamical equilibrium 
morphologies of the central dust features point towards external influences being responsible for visible dust features in early type 
shell galaxies. Inner shells are able to manifest themselves in relatively old shell systems. 

Key words. Galaxies: elliptical and lenticular, cD; Galaxies: photometry; Galaxies: interactions 



1. Introduction 



Shell galaxies (Malin & Carter 1980) have long been recognised 
as useful laboratories for learning on both the formation pro- 
cesses and the internal structure of elliptical galaxies. Soon after 
their discovery, shells were identified as tracers of "the splat- 
ter produced by a merger" (Schweizer 1980), more specifically 
minor mergers. Hence shell galaxies provided candidate config- 
urations to investigate details of the accretion process, such as 
the nature of the accreted galaxy, the dominant types of accre- 
tion orbits, the radial distribution of the accreted matter, and the 
connection of accretion events to AGN activity. If shells trace 
specific orbit configurations of accreted stars, then they poten- 
tially contain information on the three-dimensional shape of the 
galaxian potential. 



Prieur ( 1990 1 recognised different morphological categories 



of shell galaxies. Type I shell galaxies have shells antisymmet- 



* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are asso- 
ciated with program G09399 and GQ9427 



rically (interleaved) aligned along the major axis. Type II shells 
are placed all around the galaxy. Type III shells show both or 
irregular features. Numerical simulation work in the eighties 
and early nineties provided the most widely accepted framework 
for interpreting these shell morphologies in terms of mergers. 
Quinn's phase-wrapping formalism provides an elegant expla- 
nation for Type I, interleaved shells. Phase wrapping recognises 
the discrete distribution of turning points for radially-injected 
stars that oscillate back and forth in the galaxian potential (Quinn 
1984); the shells themselves are the loci of the turning points, 
where pile-up leads to increased surface brightness; these loci 
move outward in the galaxy as a consequence of the proportion- 
ality between orbital period and apocentre distance. More com- 
plex shell systems of Types II and III can be produced by minor 
mergers from non-radial orbits, on non-spherical parent galax- 
ies, or due to internal rotation in the accreted galaxy (Hernquist 
& Quinn 1988, 1989). Shells may also result from "space wrap 



ping" ( |Dupraz & Combes 1986) in the absence of radial orbit 
turn-around, when line-of-sight integration leads to an increase 
in surface brightness for debris of a satellite accreted on a high- 
angular momentum orbit. Finally, shells may result from ma- 
jor mergers between two disk galaxies, as a result of the re- 
turn of tidal tail material (Hernquist & Spergel 1992 ), whenever 
the bulge-to-disk ratio of the parent galaxies is low (Gonzalez- 
|Garcfa & Balcells|2005| ). 
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Models for shell formation not based on mergers have been 



proposed as well. These include tidal interactions ( |Thomson & 
Wright 1990; Thomson 1991) or asymmetric local star forma 



tion ( |Loewenstein et al.|1987 ). Several observational diagnostics 
may be used to test the various theories. Shells are mostly ob- 
served in isolated environments; this may indicate either a lower 
formation rate or a shorter lifetime in denser environments; it 
may indicate younger ages for shell systems in empty environ- 
ments ( |Colbert et alJ|20T)I] l. |Forbes & Thomson| fiW2\ noted 
that almost all galaxies which contain a kinematically decou- 
pled core (KDC), also show shells, suggesting a relation between 
these galaxy properties. These two properties give support to the 
merger/accretion model for shells. 

How close to galaxy centres do shells exist is a matter of cur- 
rent interest. Because shell detection requires some form of un- 
sharp filtering, whether digital or photographic ( Malin & Carter 
1980| l, and that works best when brightness gradients are not 



pronounced, shells have more often been detected in the outer 
parts of galaxies. However, inner shells contain useful informa- 
tion on the shell-making process. The existence of inner shells 
requires that orbital energy and angular momentum be removed 
from the accreted stars before these are released into the poten- 
tial of the larger galaxy, therefore inner shells imply that dynam- 
ical friction operated, and that the accreted galaxy disrupted late 
into the accretion process. Because shells trace accreted material 
with orbit apocentres at the shell radii, inner shells would pro- 
vide strong evidence for the late build-up of the inner regions of 
ellipticals through accretion of small galaxies. 

To improve on our ability to detect inner shells, we used the 
ACS camera on board HST to image the inner parts of six well- 
known shell galaxies (see Table [TJ. The spatial sampling of the 
ACS, six times better than typical ground-based cameras, cou- 
pled to the absence of atmospheric blurring, provides for a more 
accurate modelling of the underlying galaxy brightness distribu- 
tion, and a more accurate mapping of the shell profiles them- 
selves. In this paper we present photometric data in V and / for 
the galaxies and the detected shells. 

The HST images also allow a precise determination of shell 
colours. The latter may provide useful diagnostics on the vari- 
ous shell formation models. The interaction model predicts sim- 
ilar colours for the shells as the host galaxy, whereas significant 
differences in shell colours are possible in the merger models. 
To date, observations give a confusing picture on shell colours. 
Examples are found of shells that are redder; similar, or bluer, 
than the underlying galaxy. In some cases, different authors re- 
port opposite colour differences (shell minus galaxy) for the 



same shell; we mention specific instances of this in Sect 5.3 



Colour even seems to change along some shells; examples are 
NGC 2865 dFort et al.| [19861 , NGC 474 ( jPneurl [T995j ), and 
NGC 3656 (Balcells 1997| l. Errors in shell colours are very sensi- 
tive to the correct modelling of the underlying light distribution. 
HST images allow for a detailed modelling of the galaxy light 
distribution, especially near the centres, and should provide in- 
creased accuracy in the determination of shell colours. 

Another important issue is the properties of the visible dust 
in the centres shell galaxies, which might say something about 
the dust visible in the centres early type galaxies in general. Our 
observations may help to learn more about dust origin and for- 
mation theories (Lauer et al. 2 005 1 ). 

The paper presents a photometric analysis of the shells in the 
six shell galaxies imaged in our HST/ACS program. The systems 
observed contain three type I galaxies (NGC 1344, NGC 3923 
and NGC 5982) and two type II (NGC 474, NGC 2865) and one 
type III (NGC 7626). The HST/ACS images are analysed by ap- 



plying the technique of Voronoi binning (Cappellari & Copin 
2003) on the single passband images, which yields high-S/N 
brightness and colour maps to see if the shells influence local 
colour gradients.. In an companion paper, the properties of the 
globular clusters were analysed (Sikkema et al. 2006, hereafter 
Paper I). 

The data reduction is briefly summarised in Section 2. In 
Section 3 we describe how we obtained global parameters, the 
shell fluxes and colours, the production of colour maps and dust 
properties. In Section 4, we compare our observations with N- 
body simulations of shell galaxies using different models, dis- 
cuss implications of shell colours and analyse the dust properties 
of shell galaxies. We summarise our main conclusions in the last 
Section. 



2. Observations and Data Reduction 

The six shell galaxies were observed with the ACS_WFC cam- 
era between July 2002 and January 2003 with the filters F606W 
(V-band) and F814W (I band) in CR SPLIT=2 mode . The cam- 
era contains two CCDs of 2048 x 4096 pixels, each pixel having 
a size of 0".049 pixel -1 resulting in a field of view of 202" x 
202". Exposure times were on average 1000s. The inner 24 pix- 
els of NGC 2865 and the inner 8 pixels of NGC 474 were satu- 
rated in both V and I. Table 1 contains the main characteristics 
of the galaxies like Right Ascension and Declination at Epoch 
J2000.0, the extinction coefficients, exposure times and adopted 
distances throughout this paper. Detailed information about the 
data reduction can be found in Paper I. 

In addition we found three B band (filter F435W) observa- 
tions from July 2003 of NGC 7626, associated with program 
G09427 in the HST archive with exposure times of 2620, 2620 
and 2480 seconds.. We used the standard reduced images and 
combined them to remove most cosmic rays. We used a galactic 
extinction value of 0.313 (Schlegel et al.| 1998 1 and applied the 
following transformation formulae ( Sirianni et al.|2005] l: 



Bj = m(F435W) + 25.709 + 0.108(B - V) yc - 0.068(B - y)J c (l) 

Vj = m( J F606W) + 26.410 + 0.170(fi- V). 7C +0.060(B- V)^ c (2) 

Here =m(F435W) and m(F606W) are the ACS instrumental 
magnitudes and B and V are in the Johnson-Cousins system. 



3. Data Analysis 

3. 1 . Global parameters 

Information about the morphology of the galaxies and location 
of possible shells was obtained by using the ellipse fitting task 
GALPHOT (see J0rgensen et al. 1992| l; it returns information 
such as ellipticity, position angle, surface brightness and the 
C3,C4,S3,S4 coefficients ( |Carter|1978) , all as a function of ra- 
dius. A model galaxy subtracted residual image is returned as 
well, which is shown in Figures A. 1 to F. 1 for each galaxy. 

In the GALPHOT processing, background galaxies, point- 
like objects, dust features, bright pronounced shells and addi- 
tional bad data were masked out by hand in an iterative way. 
Remaining faint shell structures, having a brightness typically 
not more than 5% of the underlying galaxy emission, do not no- 
tably affect the results. The best fits were obtained by allowing 
the centre, position angle and ellipticity as free parameters to 
vary. In two cases: NGC 2865 and NGC 5982 the central regions 
do not have reliable fits (see Figures C.2 and E.2 respectively). 
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Table 1. Properties of six shell galaxies. Columns (1-11) give data from the literature (with columns (1-10) from paper I and column 
(11): de Vaucouleurs et al.| 199 columns (12-14) present data derived in this paper: (1): Galaxy name, (2): shell galaxy type, (3- 



4): Right Ascension and Declination in degrees, (5-6): Galactic extinction in V and I (in magnitudes), (7-9): Morphological type, 
apparent V band magnitude and absolute V band magnitude. (10): velocity dispersion cr in km/s, (11): effective radius in arcsec. 
(12-13): Calculated background values of the ACS V and I images in counts. (14-15): fitted Sersic index n, starting from radius r 
(arcsec). *: no stable fit possible. Fits are drawn in Figures A.3-F.3. 



These regions correspond to rough circles with diameters 24" 
for NGC 5982 and 4" for NGC 2865 respectively. The pixels 
within 0.3" of the centre of the latter galaxy are saturated. 

Global surface brightness profiles were obtained by plotting 
the surface brightness for each fitted ellipse as function of ra- 
dius. The outer parts of the profiles are severely influenced by 
the background uncertainty. It is difficult to determine reliable 
background values from the ACS images themselves, since the 
galaxies fill the whole field. Fortunately, for all galaxies, apart 
from NGC 5982, we found optical wide field data in the R band 
in the ESO archive from the Wide Field Imager (WFI) at the 
2.2m ESO/MPI telescope. The WFI camera has a field of view 
of 34' x 33', much larger than the galaxies. 

We used the ASTRO- WISE system [^( fValentijn & Kuijken 
2004) to reduce the WFI images. After subtracting a constant 
background value from the WFI images, GALPHOT was ap- 
plied on them. The ACS background values were determined 
by matching ACS surface brightness profile to the WFI sur- 
face brightness profiles. For the galaxy without WFI data, NGC 
5982, we assumed that the surface brightness profile follows a 
straight line from a certain point in a r'^ 4 -log(I) plot. The calcu- 
lated values for the background are listed in Table 1 . The final 
GALPHOT results describing the morphology of the six galaxies 
are shown in Figures A3-F3 in Appendix A to F. These Figures 
also include the results of the isophotal analysis on PA. and el- 
lipticity from WFI ground based images. The errors in V-I (top- 
right panels) were determined using the scatter in the WFI back- 
grounds. We fitted Sersic profiles to the I band surface brightness 
data using of the equation: 



yu(r) = n e + c„ 



- 1 



Such 'bulge-disk' decomposition of the profile led Schombert & 
Wallin] ( |1987| l to argue that NGC 474 is a face-on SO, a point 
which we address in Sect. 14.11 



Residual images are obtained after subtracting the galaxy 
models obtained by GALPHOT and are shown in Figures Al- 
Fl. They were solely used to identify and locate the shell fea- 
tures, but not to determine the brightness of the shells: the resid- 
ual images still show some large scale fluctuations in the back- 
ground which will disturb measurements of faint shell fluxes sig- 
nificantly. This makes it difficult to obtain reliable shell bright- 
ness from these images. A better approach, described in the next 
Section, is to work locally, within wedges. 

Except for NGC 1 344, global isophotal analysis of our galax- 
ies with other data has been done before: NGC 474 (ground- 
based B and V: Pierfederici & Rampazzo 2004; ground based 
B and R : |Turnbull et al.| | 



1 1999) 1, NGC "2865 (ground based B, 



(3) 



VandI:|Reid et al.||1994| >, NGC 3923 (ground based B and R: 
Jedrzejewski||1987]i; NGC 3923, NG C 5982 and 7626 (ground 
based V, R and I: |Bender et al.|1988) and NGC 5 982 and NGC 
7626 (HST WFPC2, V and I: |Carolk> et al.|1997[ NICMOS 1.6 
fim: |Quillen et al. 2000). As mentioned before we can also com- 
pare with WFI archive data for NGC 474, 1344, 2865, 3923 and 
7626. Comparison with these data give similar results. 



3.2. Shell radii 

Shell positions were determined by two of us (DC and GS) by 
visual inspection on an image display of the residual images de- 
scribed in § 3.1 Shell positions are listed in Table 2. Following 
Prieur, we list radii corresponding to the outermost edge of each 
shell. 

Shell radii are discussed in § |5.1| 



with c n = 2 5 (0.868rc - 0.142) (valid for 0.5 < n < 16.5; 
Caon et al.|1993 1. The fit was made from that particular radius, 
chosen by eye, where the, sometimes visible, inner plateau, will 
not disturb the fit. These starting radii are drawn in Figures A-F.3 
as vertical dashed lines. For two galaxies, NGC 474 and NGC 
2865, no stable fits were possible, with large variations for n 
depending on the starting point fitting radius r. The fitted values 
for n and starting points are listed in column (13-14) of Table 1. 

For NGC 474, adding an outer exponential to the fitting func- 
tion significantly improves the surface brightness profile fit, re- 
sulting in a smaller and less concentrated spheroid than that 
listed in Table [T] r e = 6.6" and n Ser = 3.0, and B/D = 0.71. 



www.astro-wise.org/portal 



3.3. Shell fluxes 

We have developed a three-step procedure to determine the shell 
brightness. Common to all the steps is that we work locally, 
within wedges, that are carefully placed over parts of the shells. 
For a particular galaxy the same wedges, shown in Figures A. 1- 
F. 1, were used in the two passbands. Next follows a detailed de- 
scription of the procedure. 

1. Determining the local surface brightness profile of the 
galaxy. 

First the surface brightness profile within the wedge was de- 
rived. Data points of the curve were calculated by averaging 
the pixel values within the partial elliptical rings covered by 
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Table 2. All shells identified by eye in our GALPHOT residual images. (1) Shell number/label. Asterisks indicate shells which have 
a colour determination (see also Table 3 and the appendix); (2) Shell direction relative to the centre; (3) Semi-major axis of outer 
shell border in arcseconds as measured from the galaxy centre; (4) Similar to column 3 but now in terms of effective radius; (5) 
Comments: additional information. The comments for NGC 3923 give Prieur (1988) labels + positions. 



the wedge. The elliptical ring segment had a thickness of 2.5 
pixels, with a fixed central point and an ellipticity and posi- 
tion angle, whose values were taken by averaging the I-band 
GALPHOT results of the outer galaxy regions. Pixels within 
the wedge belonging to GALPHOT masks (see Section 3.1) 
were not used. Remaining pixel outliers in the elliptical ring 
segment were removed by iterating 10 times over the set of 
pixel values, each time applying a 4 <x clipping method. The 
resulting surface brightness profile valid for the wedge is 
used in the next step. 
2. Making a local galaxy model. 

Several smooth curves were fitted to the surface brightness 
profile using Legendre polynomials of different degrees, ap- 
plying the IRAF tool 'CURFIT'. Usually the shells are al- 
ready visible as small bumps in the profile. Existence of 
shells is double checked by inspecting the residual images 
of GALPHOT. Another check is made by inspecting if the 
bumps are visible at the same position in both passbands. 



The data points in the profile which are part of the bumps in 
V and I, were left out in the fitting procedure. The fitted pro- 
file gives the model flux values in the partial elliptical rings 
mentioned in the previous step, which enables us to construct 
a 2d model image, which is valid locally for the wedge. 
3. Obtaining shell surface brightness. 

After subtracting each model image from the galaxy image, 
the residual images show the shells in each wedge for each 
Legendre model. The same recipe, as described in the first 
step, was applied again to obtain surface brightness profiles 
for each residual image, but now using the ellipticities of the 
shells, which, in case of type I shell galaxies (NGC 1344, 
NGC 3923 and NGC 5982), are much rounder than the un- 
derlying galaxy (this was already known, see e.g. Prieur, 
1988). The surface brightness profiles now clearly show the 
shell fluxes and the background is close to zero. The results 
are shown in the top panels of the Figures (i.e. 4, 5, 6 etc.) 
in Appendix A-F. Reliable shell fluxes could not be obtained 
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for left-side (North-East) of NGC 3923 and NGC 5982, be- 
cause they are not well defined features, which is related to 
their low contrast and S/N w.r.t. the galaxy. 



3.4. Shell colours 

By combining the results of two passbands for each Legendre fit, 
we obtain the average colours in the shell regions. Due to the dif- 
ferent fits (i.e. each fit is constructed using a particular Legendre 
polynomial), the derived colours show variations which increase 
for fainter shells. Consequently, we use only those shell regions 
for which all fits showed a stable answer. These regions are indi- 
cated by vertical lines in the top panels of Figures 4, 5, 6 etc. in 
Appendix A-F and typically contain data points with at least 10 
counts. The final shell colour within a region is calculated by ap- 
plying a weighted average, using the values for the colours and 
their errors derived for each Legendre fit. The resulting values 
are plotted in the Figures and listed in Table 3 . 



3.5. Galaxy colours 

Global colour maps of the galaxies were obtained by using the 
adaptive binning algorithm (Capp eTIari & Copin|2003 1. This al- 
gorithm bins two-dimensional data to a constant signal-to-noise 
ratio per bin by calculating a Voronoi tessellation. For objects 
with large gradients in S/N, as is the case for galaxies, this will 
result in smooth 2D colour images in regions with low S/N and 
therefore will show the colour gradients better than in the tra- 
ditional way. A practical example and explanation of this algo- 
rithm is given in Ferreras et al. ( 2005| l. 

The procedure to obtain the colour maps consists of three 
steps: 

First, the adaptive binning algorithm is applied on the 3x3 
binned non-background subtracted I band images by using a S/N 
of 250 and leaving out masked regions. The resulting Voronoi 
tessellation is also used for the V band image in the next steps. 
The Voronoi tessellation is further processed by applying a 
Delaunay triangulation to the central positions of the Voronoi 
cells, assigning the average flux values within the cells to these 
central positions. Finally, V-I colour maps are obtained by sub- 
tracting the appropriate background values and combining the 
images by applying the appropriate transformation formulae (see 
Paper I). 

In the outer regions, where the uncertainties in the back- 
ground become important, we still see large variations. These 
outer regions also show discrete offsets between different quad- 
rants which increase when going outwards. This is due to ran- 
dom variations of the subtracted bias level as measured in the 
overscan versus the actual bias level in the science image^ 
These differences can be as large as a few counts, which wil 
show up especially in the V-I colour maps with low signal: the 
outer galaxy regions in our images. Unfortunately, appropriate 
calibration data to correct for this effect only exists for ACS ob- 
servations later than November 2004, much later than our obser- 
vations. 

NGC 474 and NGC 3923 show elliptical red rings which 
are approximately 0.05 mag. higher in V-I than its surroundings 
(Figures 1 and 2 respectively). We believe that these are caused 
by artifacts, probably reflections, within the optical system. We 
double checked this for NGC 3923 by constructing a V-I colour 
map from ground based VLT-FORS2 images in Bessel V and I. 
The colour map does not show the ring which is illustrated in 




http://www.stsci.edu/hst/acs/documents/isrs/isr0407.pdf 



Fig. 1. Voronoi binned colour image of NGC 474. The ring is 
probably an artifact; offsets between quadrants are also visible. 



the lower panel of Figure D.4 showing, as a green line, the V-I 
profile of the ground based data: the bump due the ring between 
r=40" and 60" is not visible. 

Any correlation between shells and integrated V-I colours is 
checked by calculating the Voronoi colours in the same wedges 
as used in the previous section; in general the shell fluxes are 
so low that little of it is reflected in the integrated V-I colour 
profiles. This can be seen in the lower parts in Figures 4, 5, 6 
etc. in the Appendix A-F. We also see that most shell V-I colours 
are usually similar or sometimes redder than the colour of the 
galaxy. The signature of the rings of NGC 474 and NGC 3923 is 
also visible in these Figures as shallow, large scale, bumps with 
an amplitude about 0.05 at 25 arcsec and 50 arcsec for NGC 
474 and NGC 3923 respectively. For NGC 7626, where we also 
have B data, we only could determine a reliable shell colour for 
shell 2, which has (B - V) gahxy = 0.99 + 0.02 and (B - V) sheU = 
1.11 ±0.05 (see Figure 6). 

3.6. Shell radial profiles 

Figure 3 shows profiles for 21 bright shells in several galax- 
ies. Some shells show plateaus and are asymmetric: they 
reach a maximum flux near the outer shell border, often 
accompanied with a sudden sharp drop (examples: N474_6, 
N474.13, N2865^, N3923_8, N3923J0, N7626_2, N7626_3). 
Other shells have a symmetrical Gaussian-like shape (N474_7, 
N474.il, N3923J, N3923_15, N1344J3, N1344.6). Some 
shells seem to have double peaks (not shown in Figure but see 
the Appendix). A quantitative characterisation is obtained by ap- 
plying a Gauss Hermite fitting procedure (van der Marel & Franx 
1993) using five free parameters: y, Rq, <t, H3 and H4. The fit- 
ted values are given in Table 3 and the results are drawn as red 
curves in Figure 3. 

The shells we see are the result of a projected 3D density 
distribution in our line of sight. To get more information about 
the real or intrinsic three dimensional density distribution, we 
integrated along the line of sight using two different simple 3D 
density shell models assuming spherical symmetry and assum- 
ing a opening angle <p in our line of sight. Two other parameters 
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Fig. 3. Shell profiles with different fits. Horizontal scale in arcsec. Vertical scale in counts. The red curve is a Gauss Hermite fit to 
the data giving a quantitative characterisation of the shell shapes. Green and blue curves are the results by integrating along our line 
of sight a spherical symmetrical density model using a (ro - r) and Gaussian distribution respectively. 



used in both models are the points ro where shells have their 
maximum stellar density po. 

The first model assumes that the shell has an intrinsic 
(ro - rT 05 density distribution which is predicted by theory 
to describe the inner parts of phase wrapped shells (Dup raz &| 
Combes 1986, Prieur 19881. The second model is supposed to 
describe spatial wrapped shells ( Turnbull & Collett|1999[ ). It as- 
sumes that shells are have an intrinsic Gaussian shape with thick- 
ness r g . However this model lacks a physical basis. The results 
are again given in Table 3 and shown in Figure 3: here the green 
and blue curves depict the results of the first and second model 
respectively. In most cases the second model describes the data 
best. The 'noise' in the green curve is due to the fact that a dis- 
continuous model (the model has a sudden drop in flux to zero at 



ro) is fitted to the data: the noise in the model reflects the noise 
in the data. 



4. Notes on the individual systems 

4.1. NGC474 

The low v/<T ratio (this defines the amount of rotation), of 0.18 
(Rampazzo et al. 2006) is consistent with NGC 474 being ei- 
ther a near face-on disk, or a near spherical galaxy. NGC 474 
is part of a small evolving poor group (Rampaz zo et al.|2 006 ). 
Its shell system is very complicated as shown in Figure 1 of the 
multi- wavelength study of Rampazzo et al.. In addition to the 
shells, a loop is visible, which heads east-west, starting from a 
comma shaped feature (outside the ACS field of view), passing 
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Fig. 2. Similar as Figure 1 but now for NGC 3923. 



the sideline towards the galaxy centre. In HI, there are signs of 
tidal interactions with the nearby (at 5.5') regular spiral galaxy 
NGC 470. The same study does not make clear, however, if NGC 
470 is responsible for or even related to the shell system. 

Being classified as a type II shell galaxy, this system has been 
used to test predictions of the weak interaction model (WIM; 
Tho mson & Wright|1990 Thomson] 1 99 Here, shells are in- 



duced in the outer parts of the host galaxy due to tidal effects 
resulting from a fly-by of another galaxy. Morphologically, the 
WIM simulations show the shells as almost complete windings 
or spirals around the centre, when looking face-on. Looking at 
our GALPHOT residual image of NGC 474, however, the shells 
look more like short arcs. The shell positions and shapes (see the 
residual image Figure A.l) resemble more those of the results 



of the merger simulations shown in e.g. Figur e 5 of [Dupraz & 

!9J 



Combes ( 1986 ) , model 4 ofp ernquist & Quinn| ( fl989| l, or model 
7 of |Hernquist & Quinn ( 1988| l: all mergers of low-mass com- 
panions on non-radial orbits with a spherical or mildly oblate 
primary. In these merger simulations the shells can be very old 
(>5Gyr) and are spatially wrapped around the galaxy. 

Except for one shell, the colours of the analysed shells (see 
Table 3 and Figures A.4-A.8) are similar to the galaxy, consis- 
tent with previous studies (in B-V, V-r, R-I (Schombert & Wallin, 
1987), inB-R (T urnbull et al, 1999) and in B-V ( |Pierfederici &] 
Rampazzo 2004)). Some large shells are overlapped by differ- 
ent wedges, enabling us to compare their colours independently: 
shell 3 is located at * 40" from the centre. We splitted the shell 
into shells 3a and 3b (see Figure A.l); comparing their colours 
gives similar V-I colours: 1.21 +0.13 for 3a and 1.17±0.05 for 
3b. Similarly, shell 7 was splitted into shells 7a and 7b (located 
at about 64" from the centre, see Figure A.l). They have V-I 
colours of 0.96 +0.05 and 1.01 ±0.08 respectively. For shell 13a 
and 13b, at 103" from the centre, we find significantly differ- 
ent colours: 0.93+0.04 and 1.14+0.05. These are probably erro- 
neous values due to low shell fluxes compared to the underlying 
galaxy. Wilkinson et al. ( 2000| l also find offsets in colour up to 
0.30 mag between shell segments a similar type II shell galaxy 
0422-476. Another colour determination for this shell is given 
by Schombert & Wallin| ( |1987| ), who find redder R-I colours than 
the galaxy (galaxy R-I=0.88, shell R-I=1.09) and Pierfederici & 



Rampazzo (2004) who find slightly redder colours in B-V. The 
innermost sharp edged shell is detected at about 30" from the 
centre in the South direction. 

Of all shells analysed in this work, the only really blue shell 
relative to the integrated galaxy colour is found in this galaxy, 
which is shell 5. Interestingly, (Turnb ull et al.|1999| , found that 
their only really blue shell w.r.t. the integrated galaxy colour is 
the comma shaped feature at the SW (beyond the field of our 

We note 



ACS images, but see Figure 9 Turnbull et al. 



1999! 



that the position of our blue shell 5 lies exactly on the tail or 
loop, which is connected to the comma shaped feature. We there- 
fore suggest that shell 5 is related to this feature. Blue shell 
colours in other shell galaxies have been found in young, gas- 



rich merger remnants such as NGC 3656 (Balcells 1997 1 and 
Arp 230 (McGaugh & Bothun, 1990) as well as blue tails in 



many other interactions (see e.g. Schombert et al. 1990). It is 
therefore tempting to conclude that this entire feature is the re- 
mains of a recent small merger unrelated to the rest of the red 
shell system. 

Isophotal analysis of NGC 474 (Figure A. 3) shows that the 
ellipticity is changing fast from 0.08 to 0.24 between 10" and 
20". and back to 0.08 beyond 20". At the same radii, the posi- 
tion angle is changing from to about 20 and back. The galaxy 
contains several pronounced dust lanes within the inner 15", not 



seen previously (see, e.g., 


Ravindranath et al.|2001||Sarzi et al. 


2006). 


Ravindranath et al. 


find a point source in the centre. The 


top rig 


it panel of Figure A.3 shows that this source is 0.05 mag. 



bluer in colour than its surroundings (V* — / = 1.41). NGC 474 
contains the largest visible dust mass of our sample 1O 4 M , 



Table 4) and its centre shows peculiar kinematic behaviour ( Hau 
|etal|1996D . 

4.2. NGC 1344 

The shells in this galaxy are supposed to be the result of phase 
wrapping, since NGC 1344 is a type I shell galaxy. The colour 
of one outer shell of this type I shell galaxy was determined by 



Carter et al. 



1983 1. This shell appears to be somewhat bluer than 
the main body of the galaxy. In our GALPHOT residual data the 
innermost shell is visible at about 27" North from the centre. For 
the first time, we determined positions and the colours of some 
inner shells for this galaxy. The shell positions (Table 2) seem 
to show an interleaving pattern, although not as regularly as in 
NGC 3923. The V-I colours of the Northern shells are redder 
than the galaxy while they are similar or slightly bluer than the 
galaxy on the South side. The redder colours on the North side 
are probably due to dust, since these shells are much nearer to- 
wards the centre where dust is more present. The blue colour of 
one of the Southern shells is at least consistent with the earlier 
finding of Carter et al. The low number of shells and blue colour 
of some of these may be evidence for a relatively recent merger 
event. 



4.3. NGC 2865 

This galaxy is classified as a type II shell galaxy. However, look- 
ing at our residual images this system looks more as if it belongs 
to the type III class, with lots of irregular features and loops. 
The core is much bluer (V-I « 1.05) than the outer parts (V-I « 
1.15). This blue colour is related a young stellar population (0.4- 
1.7 Gyr) which forms a KDC ( |Hau et al.|1999). NGC 2865 also 
contains an incomplete HI disk (Schimino vich et al.|1995 1. An 
interesting result is that a bright HI patch is coinciding with a 
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a [arcsec] 

Fig. 4. Light profile of the galaxy along the wedge containing the 
brightest shell 2 of NGC 2865. The shell is about 0.8 magnitude 
brighter in V than the galaxy. 



bright shell observed earlier (labelled as shell 2 in this work and 
called shell 2B by Fort et al. ( 1986)). Figure lc of Schiminovich 
et al. shows an overlay of the HI data and Fort's schematic shell 
map. If these are related, this means that the shell is moving to- 
wards us, implying a spatially wrapped shell and confirming that 
this is not a type I shell galaxy. However, the asymmetric shell 
profile seems more consistent with a phase wrapped shell. Fort et 
al. derived a slightly redder colour for shell 2 of ( V - R) j hns<m= 
0.84 +0.09, compared to the galaxy colour (V-R)j=0J5 +0.03. 
A similar colour difference shell 2 and the galaxy is found in this 
work: V-I=l.ll ±0.05 vs local galaxy V-I=0.98 +0.01. The 'jet', 
as reported by Fort et al. looks more like a loop (see top of Figure 
C.l), which may or may not be connected with the bright shell 
and HI. 

Compared with the simulations, the residual image looks 
similar to situations in for instance Figure 2 or 6 of Hernquist 
|& Quinn| ( fl988]l, Figure 10 of |Hernquist & Quinn| ( | 1989) 1, or 
Figure 11 of |Dupraz & Combes ( 1986] >. All these simulations 
use small disk galaxies as the intruders, which is supported by 
the presence of an HI disk in NGC 2865. 

The bright shell 2, coinciding with the HI patch, is probably 
the best candidate for follow-up spectroscopy of our whole sam- 
ple, because of its high contrast with respect to the galaxy light 
(a;0.8 mag. higher in V than the galaxy, see Figure 4). 



4.4. NGC 3923 

This is probably the most studied shell galaxy and the proto- 
type of the Type I shell class, where shells are placed in inter- 
leaving order along the major axis. Prieur ( 1988), using ground 



we list and compare his shell positions with our results. There is 
a good agreement; small offsets between positions are likely due 
to measurement errors. Fort et al. (1986) give colours for three 
shells in NGC 3923, with only one of them in our field of view 
(our shell 8): his shell 1 has a colour of {V-R) Johnson =0.68 ±0.18 
with the local galaxy colour being about ( V - R)j = 0.82 + 0.03. 
Our results for this shell are V-I=1.29 +0.05 and local galaxy 
colour of V-I=1.28; the colour difference being consistent with 
Fort' s work. On the West side, we find that all shells have a sim- 
ilar or slightly redder colour than the galaxy. We were not able to 
determine reliable shell colours on the East side. The local mod- 
els, using different degrees for the Legendre polynomials, do not 
give stable answers. This is related to the low S/N of the images 
and badly defined shell features. 



4.5. NGC 5982 

NGC 5982 is a Type I shell galaxy (Figure 5). The galaxy is well 
known for its KDC (Wa gner et al.|1988] ), recently confirmed us- 
ing 2D central mapping using OASIS ( |McDermid et al. 2006 



and SAURON (Em sellem et al.|2004| l. Analogy with the KDCs 
of Emsellem et al. (2004) indicates that the KDC is probably 
a rotating central feature, i.e. a disk. The ellipticity becomes 
very round in the inner 2" of the galaxy. We were not able to 
smoothly fit the central regions. Even after using 4 harmonics, a 
quadrupole with wings containing a flux of about 200 counts (+ 
1-2 % of the galaxy flux), is still visible. The C4 coefficient in 
this region is about -0.02 (Figure E.3) indicating a boxy structure 
(Carter|1978|l. 2D spectral mapping of this region (Emsellem et 



based CCD data and photographic plates, mapped the whole 
shell system. A comparison shows that we do not find any other 
shells than those already given by Prieur, despite our much better 
resolution and galaxy subtraction near the centre (compare our 
Figure D.l with his Figures 1-5). In the last column of Table 2, 



|al.|2004| l show that the stellar velocities exhibit a 90 degree off- 
set in the central region with respect to the outer regions. 

The bad GALPHOT fit limits us in giving a final answer 
about how far the shells do extend to the centre. The innermost 
detectable shell is located about 8" East. Assuming the galaxy is 
about two times further away than NGC 3923 (see Table 1), the 
distance from the centre of this shell is comparable to the inner 
shell of NGC 3923. Double checking using other methods (un- 
sharp masking and GALFIT (Peng 2004), which fits symmetric 
2D models), show no more inner shells. 

Unfortunately reliable shell brightness could not be deter- 
mined, because of their faintness and low contrast. A deeper, 
3600s, ground based image, was found in the ING Archive 
(taken with the 2.5 meter Isaac Newton Telescope INT in July 
1989 in the R band). After again using GALPHOT and subtract- 
ing galaxy-model, we obtain Figure 5, which shows the shells 
with somewhat higher contrast. The outer shell 24 on the minor 
axis is barely visible in the ACS images and is 90 degrees dis- 
placed with respect to the inner shells . The next bright outer 
shells 20 and 21 are slightly displaced with respect to the inner 
shells. This peculiar shell morphology resembles the scenario 
shown in Figure 6c of Dupraz & Combes ( 1986). In this simu- 
lation a small elliptical falls into the potential of a prolate E3.5 
galaxy with an impact angle of 90 degrees. The same misalign- 
ments occur for the outer shells in this event. 

This galaxy is classified as an YP galaxy (i.e. containing a 
young stellar population similar to NGC 2865) with a fine struc- 
ture index £ 2 of 6.8 ( |Michard & Prugniel||2004| l. Michard & 
Prugniel do not mention shells. Looking at our images we es- 
timate that there are at least eight shells, which would imply a 
much higher £2 of about 1 1 . 
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Fig. 5. Ground based residual image of the shell system NGC 
5982, with shell labels (see Table 2). This morphology resembles 
the scenario shown in Figure 6c of Dupraz & Combes (1986 ), 
where a small galaxy has fallen in the host galaxy at an offset 
angle of 90 degrees. 



4.6. NGC 7626 



Its core kinematics (Bal cells & Carter|1993[ ) and the fact that it 
has bright globular clusters (Paper I) indicate a (minor) merger 
event which happened one or several Gyr ago. This is the first 
time that shell colours and brightness have been obtained for this 
galaxy. A shell on the East side was detected by Jedrze jewski| 
|& Schechter| (1988 ); this shell lies outside our field of view. 
Another possible shell lying on the S.W. side, was detected by 
|Forbes & Thomson| ( |1992[ ). This shell is also seen in our data. 
We only detect one other, very sharp, edged shell to the N.E. The 
structure of these shells looks somewhat like those of the simu- 
lations shown in Figure 5 of Dupraz & Combes ( 1986 ). Here a a 
small spiral (1% mass of host) was thrown into the potential of 
an E3.5 oblate galaxy. Frame 3 (after 4 Gyr) of this Figure looks 
very similar to the NGC 7626 shell system. The zero rotational 
momentum encounter should create phase wrapped shells which 
is probably evidenced by the radial shell shapes (see Section 
5.2). Further support for a disk intruder galaxy comes from the 
shell colours. In both (B-V) and (V-I), the brightest, inner shell 
is redder than the galaxy (see Figure 6), which is probably due 
to dust (see Section 5.3). The fainter outer shell has a similar or 
slightly bluer (V-I) colour than the galaxy. 



5. General results 

5.1. Shell radial distributions 

As discussed in Sect. [T] inner shells contain useful clues to the 
shell formation process, and indeed, one of the central goals of 
the HST-based imaging program was to determine how close to 




Fig. 6. Colour-colour, B-V vs V-I, diagram of the brightest, in- 
ner, shell in NGC 7626. Cross + errorbars represent the local 
galaxy colour, while the triangle with errorbars represents the 
derived shell colour 



the galaxy centres we find shells. Our determined shells and shell 
radii (see Sect. |3.2) > are listed in Table 2, given in arcsec as well 
as in units of the effective radius r e (from col. 14 of Table 1). 

Innermost shell radii span a wide radial range, from r m ; n /r e = 
0.24 in NGC 5982 to r min /r e = 2.86 in NGC 2865. We find a 
slight tendency for type-I shell systems to extend within r e and 
for types II and III to lie in the outer parts: the two innermost 
shells, in units of r e , are in NGC 5982 and 3923, two Type- 
I galaxies, while the third galaxy with r m \ n /r e < 1, NGC 474 
(type-II), has an uncertain value of r e : as shown in Sect. |3.1| r e 
might be significantly smaller than listed in Table 1; adopting 
such smaller value for r e , we would get r m ; n /r e = 4.2. However, 
not all type-I shell systems show inner shells: for NGC 1344 we 
find r min /r e > 2. 

The HS T imaging has revealed shells in the inner two kpc 
for two of the type I shell galaxies: NGC 3923: r m ; n = 8" = 
1.8 kpc (at an assumed distance of 20.0 Mpc) and NGC 5982: 
r m j n = 8" = 1.7 kpc (at an assumed distance of 41.9 Mpc). 
This is interesting as it clarifies that shells do indeed form near 
the galaxy nuclei. This region is easily available to spectro- 
scopic kinematic measurements, hence the correspondence be- 
tween kinematic features and shells may be revealed. 

For the galaxies with no inner shell detections, the most 
straightforward interpretation is that shells never formed at those 
radii. This conclusion may be too simplistic. Three conditions 
need to concur for the detection of inner shells: (z) shells need 
to form; (zz) shells need to survive until the observation epoch; 
and, (Hi) they need to be detectable through the shell-detection 
methods employed. 

The second and third conditions listed above make shell de- 
tection harder as we look closer to the galaxy centres. We first 
address shell survival. We expect shell lifetimes to be shorter 
near the centres: phase mixing scales with dynamical time, hence 
it is faster near the centre; shells should lose contrast and fade 
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away faster near the centres. Furthermore, galaxy centres are 
more dynamically active than the outer parts, as any object that 
merges after the formation of the shell system and reaches the 
centre, will gravitationally perturb the orbits of the stars that de- 
fine the shells. In this respect, it is quite suprising that we see 
inner shells in NGC 3923 and NGC 5982, wh ich have a lar ge 
number of shells and are relatively old systems (Nulse nfl989| ) 

Shell detectability becomes progressively more difficult as 
we approach galaxy centres. As mentioned above, the pro- 
nounced brightness gradients in the inner regions of ellipticals 
lead to the break-down of unsharp-mask techniques. Because 
shells near the centres are closer to each other than further out, it 
may become difficult for the detection algorithm to pick up the 
underlying, non-shell brightness levels. Finally, elliptical nuclei 
are known to be dusty ( Lauer et aT||1995 Phillips et al.||1996~ 



|Peletier et aLl|1999[ |Ravindranath et al.||2001) . Dust may act 
in two ways. It may simply hide the shells: examples of such 
an effect are NGC 474, NGC 1344 and NGC 7626. Dust may 
also perturb the general light distribution, so that the smooth 
galaxy model one generates and subtracts to reveal the shells 
has too strong residuals for the shells to appear. Typically, the 
underlying light distribution needs to be smooth to within a 
few percent for the shells to show up. Strong dust patches eas- 
ily lead to stronger third- or fourth-order Fourier residuals in 
the isophotes. Examples of this situation are NGC 2865 and 
NGC 5982. Clearly, HS T imaging at NIR wavelengths would 
allow us to see through the dust and would strongly improve the 
chances of detecting inner shells in ellipticals. 



5.2. Shell brightness profiles 

Looking at Figure 3, it is clear that in general the Gaussian model 
fits the data better than the (ro - r)~ 05 model. The latter model 
especially fails for lower r. However, this model is only meant 
to describe the flux behaviour very close to the shell maximum 
at projected radius r max . The (ro - r)~ 5 model also predicts a 
fast dropping flux at distances slightly larger than r max . This is 
indeed seen at the two bright shells in NGC 7626 (see Figure 
3). Here the flux drops from its maximum to zero within a small 
interval relative to the shell size. In general this seems to hap- 
pen for shells with a plateau, i.e. N474_6, N474J3, N2865-2, 
N3923_8, N7626_2, N7626_3. 

The (ro - r)" 5 model is supposed to describe phase wrapped 
shells ( |Dupraz & Combes]|1986| |Prieur|[T988| ), but is does not 
seem to fit well most of the bright shells in type I shell galax- 
ies NGC 1344 and NGC 3923, where shells are expected to be 
the result of phase wrapping. For NGC 3923 this may be related 
to the age of the shell system. The large number of shells indi- 
cate an old age. The shells may smooth out as a result of their 
older age and will not have razor sharp edges as the (ro - r)~° 5 
model assumes. On the other hand: NGC 1344, showing only 
a few shells, indicating a younger age, does not do much bet- 
ter. NGC7626, which is probably very young (see Section 4.6), 
might be a better example where the model works. We conclude 
that the (r - r)~ 5 model works best for younger shells. Older 
phase wrapping shells probably have a more extended structure 
and density profile than the theoretical (ro - rT° 5 model, for in- 
stance due to internal velocity dispersions in the intruder galaxy. 
It seems that the (ro - r)~ - 5 is just too simplistic to describe the 
real shells. 



5.3. Shell colours 

In our galaxy sample we find only one shell with blue colours. 
All other shells have similar or redder colours. Red shell colours 
are also found by many others like in NGC 7600 (Turnbul l et al.| 
[19991 ), IC 1459 ([Forbes et al.|1995a) NGC 7010, NGC 7585 and 
IC 1575 (Pierfederici & Rampazzo 2004). The redder colours 
could be explained by at least four scenarios: 

- The stars in the shells are on average older than those in the 
main body of the galaxy. This may occur if the shells be- 
long to older parts of the intruder galaxy, e.g. the bulge, or 
if the interaction or merger led to the formation of young 
stars from gas throughout the galaxy. N-body simulations 
show that the best reproductions of shell morphologies are 
obtained by using very small intruder galaxies, with only a 
few percent of the mass of the host. It is therefore not very 
likely that they will form sufficiently many new stars which 
could lower the average colour of the whole galaxy. 

- The stars in the shells are more metal rich than those in the 
main body of the galaxy. This is highly unlikely given the ex- 
pected small mass of the intruder galaxy and usually positive 



correlation between metallicity and galaxy mass (Sandage & 
| Visvanathan| 1 978) ) . 
- The stars in the shells have different, redder, populations than 
the underlying galaxy. This scenario only works for very spe- 
cific conditions. If the progenitor galaxy is a (small) late type 
(star-forming) galaxy and star formation is truncated after 
the merger event, then after some several 10 s yrs the light 
of the original stellar population will be dominated by RGB 
and AGB stars, which will redden the integrated colours of 
the shells. This reddening effect has been demonstrated by 
Maraston (2005 ). Her Figure 27 (middle left panel) shows 



an enhanced reddening in V-I after several 100 millions of 
years, mainly due to AGB stars,. However, the amplitude of 
the reddening she found is not enough to account for our red 
shell colours, which are sometimes even redder than the el- 
liptical galaxy colour. Resolved data on the shell or stream 
in M 83 also shows significant amounts of AGB and RGB 
stars ( de Jong et al.|2 007 ), but his data are not deep enough 
to calculate a reliable global colour for this stream. 
The shells contain more dust per unit stellar mass than the 
main body of the galaxy. Here, the problem is to explain why 
shells have more dust per unit stellar mass. Several possibil- 
ities can be thought of. The first possibility is related to the 
previous item: if RGB and AGB stars make up a significant 
part of the population, their large mass loss (mostly AGB 
stars) will result in more dust residing in the shells ( Athey et 
al.|2002[ ). Another possibility is based on theory: it has been 



shown that it is possible for gas or dust to remain connected 
with the shell stars after a small merger (Kojima & Noguchi 
[1997] ICharmandaris & Combes] [20^01 . Observational ev- 
idence for significant amounts of dust residing in a shell 
was found in NGC 5128 ( |Stickel et aTJ|2004"l ). HI gas in 
shells has been found in M 83, NGC 2865 and N GC 3656 
( |Schiminovich et al.|[Y995) |Balcells et aT1[200TT ). A third, 
speculative, idea explaining the presence dust in shells, is 
that dust is swept up by the shell stars as they pass through 
the potential of the galaxy. This should be tested using sim- 
ulations. In the ISO archive, we found ISOCAM ( |Kessler et 
|al.|1996| ) data for NGC 1344 and NGC 7626. These obser- 
vations, taken at wavelengths near 9.5 fim, could in princi- 
ple detect the warm dust. Although both galaxies show red 
shells, we find no evidence for enhanced emission at the shell 
regions in these two galaxies. 
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5.4. Dust in the centres of shell galaxies 

All of our galaxies show visible dust features, mainly found in 
the central parts of the galaxy (see Figures 2 in Appendix A- 
F). Following Tran et al. (2001 1 the morphologies of the visi- 



ble dust features can be divided into two groups, i.e.: 1) nuclear 
ring or disks-like structures and 2) filamentary structures and/or 
(small) dust patches. All our galaxies show at least features of 
group 2 (filaments: NGC 474, NGC 2865, NGC 3923 and NGC 
7626; small dust patches: NGC 1344, NGC 2865, NGC 3923 
and NGC 5982). Although NGC 5 982 was listed before as a 
dust-free galaxy ( |Sarzi et aL||2006| l, we see several patches in 
the residual frames (two dust patches are visible in Figure E.2 at 
about 6.5" E and N.E.). 

NGC 7626 shows both a dust lane and a nuclear ring within 
the inner arcsec (the ring was already reported by Forbes et al. 
(1995b)). A combination of these two dust features is not seen 
very often: Lauer et al. (2005) do not find any example in a sam- 
ple of 77 early type galaxies. Saturation in the core in NGC 2865 
and a bad fit of the inner regions in NGC 5982 hinders a conclu- 
sion about the presence of nuclear dust rings in these galaxies. 
The bad fit in NGC 5982 may be related to the presence of a 
KDC in the inner regions (Wagner et al. 1988| |McDermid et 



al. 2006). When we apply GALPHOT to much lower resolution 
ground based data, we see similar (bad) results (see Figure 5). 
Lauer et al. are able to make a better fit, however they do not 
detect the dust patches and shells, probably due to smoothing 
in their modelling procedure. All dust morphologies are sum- 
marised in column two of Table 4. Column three of the same 
table lists the position angles valid for the filaments and nuclear 
disk. 

'Visible dust masses' were obtained by using the method of 
van Dokkum & Franx ( 1995). They assume that the visible dust 
acts as a foreground screen w.r.t. the background galaxy light, 
which will result in a lower limit for the dust mass and will have 
large errors of the order of 50%. The dust mass is derived with 
the following expression: 



M d =< A v > zr 



(4) 



with < A v > extinction measured at a pixel, £ the surface area 
(using distances listed in Table 1) and r = 6x 10~ 6 mag kpc 2 M Q ' 
the extinction coefficient per unit mass. Only areas were selected 
where the extinction/dust is visible by eye. Ay was calculated for 
each visible dust feature by dividing the real and model images. 
Values of Ay are all lower than 1. All derived dust masses are 
listed in column 5 of Table 4. The masses are of similar mag- 
nitude to those found by other authors (van Dokkum & Franx, 
1995; |Tranetal.|2001| i. 

5.5. Dust origin 

Currently, there are at least four scenarios which explain the 
presence of dust in the centres of early type galaxies. The dust 
survival time, which depends on the main destruction mecha- 
nism (sputtering by the hot X-ray gas), is expected to be rela- 
tively low in the centres of early type galaxies (10 7 - 10 8 yrs 
depending on the electron density, ( |Draine & Salpeter||1979| 
Tielens et al. 1994) >. A problem with these timescale calculations 
is that it does not take into account the effect of self-shielding 
in dust clouds. This may enhance the survival time consider- 
ably. The ubiquitous presence of dust in the centres of early type 
galaxies is difficult to explain without some rate of replenish- 
ment. Mathews & Brighenti (2003) showed that it is possible to 
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Table 4. Properties of dust. Column 2: Morphology of the dust 
where d=disk, f=filament, p=patchy dust; column 3 and 4: PA. 
in degrees of main dust feature and galaxy respectively; col umn 



5: dust mass in 10 M Q as determined using the method of 



|Dokkum & Franx ( 1995[ >; column 6: absolute V magnitude of 
galaxy 



form dust clouds in the centres of early type galaxies by accumu- 
lating dust from stellar winds. Other scenarios use external influ- 
ences like accretion from flybys or mergers with other galaxies. 
We will now discuss the dust properties in the shell galaxy sam- 
ple, where external influences are evident. 

Using HST archival data, about half of all elliptical galaxies 
exhibit visible dust features, which are equally present in power- 
law and core galaxies ( |Lauer et al.| 2005 47% of 177 in field 
galaxies). Assuming that a similar dust detection rate of 50% is 
representative for our galaxy sample and given the fact that dust 
is visible in all our galaxies, we can reject the statement that 
our sample belongs to the parent set of normal early type galax- 
ies with visible dust at the 97.5% level. However we should also 
take into account the possibility that our sample is biased by con- 
sidering much higher dust prevalence in certain classes of early 
type galaxies. This happens for instance in radio-loud galaxies, 
having dust detection probability of about 90% (van Dokkum & 
Franx, 1995; Verdoes Kleijn et al.|1999[). Our sam ple hosts one 
radio galaxy: NGC 7626 ( [Hibbard & Sansom|2003| l. Second, we 
consider a possible bias due to the presence of ionized gas. It is 
well known that dust is almost always accompanied with ion- 
ized gas in early type galaxies (|Macchetto et al 



|[T996) |Sarzi et 

|al.|2006| ). Due to selection effects in the detection of visible dust, 
the converse is not true, although the probability to detect dust 



in early type galaxies with ionized gas is still quite high (Tran et 
al. 2001 ). Modern instrumentation detect emission in about 75% 



of early type galaxies (Macchetto et al.|1996 Sarzi et aLf 2006 



Compared with these detection rates, our sample does not seem 
to be biased: only three galaxies (NGC 474, NGC 5982 and NGC 



7626) have low levels of Hff+[NII] emission (Verdoes Kleijn 
|et al. 1999 2002) with luminosities below the median value of 
~ r ^, determined for a nearly complete magnitude 



2 x 10^ ergs s 

limited sample of nearby galaxies ( Ho et al.]|1997 1. Combining 
these biases and dust detection probabilities still implies that we 
can reject the statement at the 95% level, and that shell galaxies 
have a higher dust prevalence than normal early type galaxies, 
contrary to an earlier finding by Sadler & Gerhard ( 1985). 

Number counts of the morphology of the dust features, oc- 
curring in our sample, also differs from normal early type galax- 
ies. Dust features in the centres of early type galaxies come in 
two types: regular rings and irregular shaped patches or lanes. 
Number counts give a ratio of 3:5 for the two types of dust fea- 
tures. They are almost never both seen in one galaxy ( |Lauer et| 
|al.|2005| l, which has been used as evidence for an (episodic) dust 
settling sequence scenario ( Tran et aL]|2001 Lauer et aT] |2005 



Verdoes Kleijn & de Zeeuw|2005| l. Here, first the irregular dust 
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patches and lanes form in some way, while some time later these 
dust features move to the centre and form a disk in dynamic 
equilibrium. All the dust features in our sample seem to be out 
of dynamical equilibrium as they show up as irregular patches 
or lanes. The probability for this to happen, assuming a regular 
to irregular ratio dust feature of 3:5, is 6%. These considerations 
lead to the conclusion that external influences are responsible for 
the ubiquitous presence of dust in shell galaxies. 

We will now briefly discuss the dust properties of the indi- 
vidual shell galaxies and see how they fit into this discussion and 
shell formation theory. 

NGC 474 contains several pronounced dust lanes within the 
inner 15" (see Figure A.2), not detected previously (see e.g. 
|Ravindranath et al.||2001| and |Sarzi et al.||2006| >. It contains the 
largest visible dust mass of our sample (« 10 4 M o ). Sarzi et al. 
report a misalignment of the stars and ionised gas by 74 ±16 de- 
grees. Comparing the position angle of the dust lanes with the 
velocity maps of Sarzi et al. show that the dust is aligned with 
the velocity field of the ionised gas and not with the stars. This 
probably implies a connection between the dust and the gas and 
is seen in many other early type galaxies (Goudfr ooij et al.[ 1994; 



Ferrari et al. 



1999; Sarzi et aL|2006| l. The velocity structure of 
the gas, the dust morphology of NGC 474 (many dust lanes re- 
siding on top of the supposed bulge) and the large offset between 
stars and gas+dust disturb the picture of this galaxy being an SO, 
and point towards an external origin for the dust. 

The well defined type I shell galaxies, NGC 1344, NGC 3923 
and NGC 5982, are expected to be created by a minor merger 
with a small, non-rotating dwarf galaxy. While the dust content 
of NGC 1344 and NGC 5982 is quite low (several time 10 2 M o , 
Table 4) and corresponds to such a scenario, the dust mass of 
NGC 3923 is an order of magnitude higher. Most (80%) of the 
visible dust in NGC 3923 belongs to a large patch, visible in the 
NE direction at 50 arcsec from the centre (see Figure D.l). This 
patch was shown to be part of NGC 3923 and also emits small 
amounts of Ha+[NII] ( Pence|[T986|l. At the projected d istance 
of this patch, the electron density (Fukazawa et al. 2006) corre- 
sponds to a dust sputtering minimum survival time of 4x 10 7 yrs. 
Another 5% resides in small patches at 4.5 arcsec from the cen- 
tre (Figure D.2). At this projected distance, the dust sputtering 
minimum survival time is « 10 6 yrs. The rest, 15%, is located 
in a long diffuse dust lane NW from centre. The amount of dust 
in this galaxy does not conform to the minor merger picture of 
a small elliptical dwarf galaxy falling into a much larger poten- 
tial (already noted by |Carter et al.|1998] l. An internal origin for 
the dust (e.g. ejection from stars) seems to be the most likely 
scenario. 

In NGC 2865 and NGC 7626, comparison of simulations 
and shell structures point towards gas rich intruder galaxies (see 
Section 4.3 and 4.6 respectively), which has resulted in recent 
(0.4-1.7 Gyr) star formation in the core of NGC 2865 (Hau et 



al. 11999)1 and likely the creation of new globular clusters (Paper 

— — ' — i „ ^ . „ i _ _ . , 



I). The dust content of several times 10 3 M Q (see Table 4) is dis- 
tributed in diffuse layers and patches near the centre (Figures C.2 
and F.2). In NGC 2865 we can again calculate the minimum dust 
sputtering survival tim e because the electron density is known 
(Fukazawa et al.|2006t. This results between 3 x 10 7 and 1 x 10 8 



years for the inner and outer dusty regions respectively. This is 
again much lower than merger timescale of several 10 s to 10 9 yr. 



6. Summary 

Using observations in V and I with the ACS on board the HST, 
we analysed the properties of shell systems, in particular their 



colours, morphologies and dust properties. For most shells listed 
in this paper, we determined their colour for the first time. For 
those shells for which their colour had already been determined, 
we find similar results, giving support to the quality of the result 
of our methodology. In general we find that colours of shells 
are similar or redder to the colours of their host galaxies. We 
attribute the red colour to dust which is physically connected to 
the shell. In some cases, a different stellar population as a result 
from a truncation of star formation, may also redden the shells. 
The only blue shell is found in NGC 474, which is very likely 
related to a long tail and probably a very recent minor merger 
event. N-body merger simulations, rather than simulations by the 
interaction model, describe best the observed morphologies of 
the shell systems. 

We detect out of dynamical equilibrium central dust features 
in all our galaxies. Comparison with a set of 'normal' elliptical 
galaxies, which have a dust detection rate of 50%, implies an 
external origin for central dust found in shell galaxies. However 
this is in contradiction with theoretical predicted dust survival 
times. Better models of dust behaviour in centres of early type 
galaxies, which include self shielding, are needed to solve this 
problem. The best shell candidate for follow-up spectroscopy 
has been found in NGC 2865. 

Innermost shells are found in the type I shell galaxies NGC 
3923 and NGC 5982 at a distance about 2 kpc from their centres. 
These shell have survived for a long time, since both galaxies 
have relatively old shell systems. 

Current models to describe the profiles of phase wrapped 
shells probably work best for young shells. 
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Fig.A.l. ACS residual image of GALPHOT for NGC 474 in 
V, surface brightness of shells were determined in the wedge 
regions using ellipticities as indicated by the white strips. The 
size of the ACS field of view is 202 x 202 arcseconds. 




Fig. A.2. Inner region residuals of NGC 474 in V (40x40 arcsec). 
Dust is directed in the NS direction and aligned with the ionized 
gas (to be compared with the Figures in Sarzi et al. ( 2006| l) 



Appendix A: Results for NGC 474 
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Fig. A.3. Morphological data NGC 474, blue and red curves rep- 
resent V and I band data respectively. Top left: surface brightness 
profiles, corrected for background (see text). The black line is a 
Sersic fit to the I band surface brightness data right of the ver- 
tical dashed line. Middle and bottom left: position angle and el- 
lipticity respectively. Green data points represent WFI data. Top 
right: Global V-I profile. Top middle and bottom: S4 and C4 re- 
spectively. 




a) [arcatc] 

Fig. A.4. Photometric data in the wedge covering shells 3b and 
9. Top panel: V band flux in counts averaged along wedge us- 
ing ellipticities as indicated by the white strips drawn into the 
wedges in Figure A. 1 . The region within the vertical dashed lines 
was used to calculate the shell V-I colours which are shown in 
the bottom panel. The bottom panel also shows the galaxy V-I 
colour profile, calculated from a Voronoi binned image. 
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G.Sikkema et al.: HST/ACS observations of shell galaxies: inner shells, shell colours and dust. 
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Fig.A.5. Photometric data in the wedge covering shells 3 a and Fig.A.7. Photometric data in the wedge covering shells 1, 7b 
6; Description: see Figure A.4 and 13b; description: see Figure A.4 
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Fig. A.6. Photometric data in the wedge covering shells 7a and Fig. A.8. Photometric data in the wedge covering shells 4, 5 and 
13a; description: see Figure A.4 11; description: see Figure A.4 
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Fig. B.l. Residual image of GALPHOT for NGC 1344 in V with 
wedges left and right. The field of view is 202 x 202 arcseconds. 

Appendix B: Results for NGC 1344 



E -1 

log(a) [arcsee] 



Fig.B.3. Morphological data NGC 1344. Description: see NGC 
474, Figure A. 3 
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G.Sikkema et al.: HST/ACS observations of shell galaxies: inner shells, shell colours and dust. 
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Fig.B.4. NGC 1344 shells 1 and 3 in the wedge on the North 
side. Description: see Figure A. 4 
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Fig.B.5. NGC 1344 shells 8 and 10 in the wedge on the South 
side. Description: see Figure A. 4 



G.Sikkema et al.: HST/ACS observations of shell galaxies: inner shells, shell colours and dust. 
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Fig. C.l. Residual image of GALPHOT for NGC 2865 in V 
with a wedge overlapping a bright shell. A loop is visible in the 



NS direction. This morphology was already drawn in Fort et al. 
( fl986| l. The field of view is 202" x 202". 

Appendix C: Results for NGC 2865. 




Fig. C.2. Inner region residuals of NGC 2865 in V (40x40 arc- 
sec). Dust is visible in the South direction. Fitting problems in 
the central region are due to saturated pixels. 
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Fig. C.3. Morphological data NGC 2865. 
474, Figure A. 3 



Description: see NGC 
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G.Sikkema et al.: HST/ACS observations of shell galaxies: inner shells, shell colours and dust. 
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Fig. C.4. NGC 2865 shell 3 in the wedge placed at the East side. 
Description: see NGC 474, Figure A.4 



G.Sikkema et al.: HST/ACS observations of shell galaxies: inner shells, shell colours and dust. 
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Fig. D.l. Residual image of GALPHOT for NGC 3923 in V with visible was also detected by |Prieur| ([1988} using ground based 
wedge. Note the large dust patch at the NE and faint dust lane data. No more other inner shells are detected. 



within the wedge in the SW direction. The field of view is 202" 



x 202". 

Appendix D: Results for NGC 3923 
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Fig.D.3. Morphological data NGC 3923. Description: see NGC 
474, Figure A. 3 
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G.Sikkema et al.: HST/ACS observations of shell galaxies: inner shells, shell colours and dust. 
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Fig.D.4. NGC 3923 shells 1, 3, 8, 10, 13, 15 and 17 in southern 
wedge. The lower panel shows shell and local galaxy V-I colours 
determined within the wedge: the points with errorbars are the 
shell colours. The solid line and green open circles represent V- 
I local colours from ACS and ground based VLT-FORS2 data 
respectively. The signature of the ring (the bump between r=40" 
and 60") is not visible in ground based data. 



G.Sikkema et al.: HST/ACS observations of shell galaxies: inner shells, shell colours and dust. 
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Fig.E.2. Inner residuals of NGC 5982 in V (60x60 arcsec). A 
Fig.E.l. Residual image of GALPHOT for NGC 5982 in V. small dust lane is yisible on the major axis in the E directiori; a 

Shells are barely visible. The field of view is 202 x 202 arcsec- more pronounce d small dust patch is visible NE. The shell labels 
onds. The shell labels are best visible in the electronic version. are best yisible in the electronic version. 



Appendix E: Results for NGC 5982 




Fig.E.3. Morphological data NGC 5982.. Description: see NGC 
474, Figure A. 3 
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G.Sikkema et al.: HST/ACS observations of shell galaxies: inner shells, shell colours and dust. 
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Fig.F.3. Morphological data NGC 7626.. Description: see NGC 
Fig.F.l. Residual image of GALPHOT for NGC 7626 in V . 474, Figure A.3 
Two bright shells are visible with the wedges overlapping them. 
The field of view is 202 x 202 arcseconds. 





Fig. F.2. Inner region residuals of NGC 7626 in V (40x40 arc- 
sec); the two shells are already outside the field of view, but the 
dust and a lot of globular clusters are clearly visible. 



Fig.F.4. NGC 7626 shell 3 in the wedge at the SW side. 
Description: see Figure A. 4 
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Fig.F.5. NGC 7626 shell 2 in the wedge at the NE side. 
Description: see Figure A.4 
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Fig. F.6. Similar to previous Figure but now the B band residual 
flux is drawn in the top panel while B-V colours are drawn in the 
bottom panel. 



